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In the hard core hmit, interacting vortices in planar type II superconductors can be modeled 
as non-interacting one dimensional fermions propagating in imaginary time. We use this analogy 
to derive analytical expressions for the probability density and imaginary current of vortex lines 
interacting with an isolated bent line defect and to understand the pinning properties of such 
systems. When there is an abrupt change of the direction of the pinning defect, we find a sinusoidal 
modulation of the vortex density in directions both parallel and perpendicular to the defect. 



I. INTRODUCTION 



The statistical physics of flux lines in high Tc superconductors has attracted considerable experimental and the- 
oretical attention. Controlling the pinning properties of the magnetic vortices that penetrate the superconducting 
material above a critical external field Hd can be the key in achieving dissipationless electrical current flow 1']. The 
interplay between point and correlated disorder, thermal fluctuations and vortex- vortex repulsion controls the flux 
line configurations and leads to a variety of different phases 0, d, [1] ■ 

Recent advances in manufacturing of high quality thin superconducting slabs and films, and experimental techniques 
for probing such systems, has made possible to observe mesoscopic vortex dynamics and image individual vortices 
confined to two dimensions Q. The trajectories of vortices in these 2-d systems can be mapped to the world 
lines of 1-d bosons and in that sense planar superconductors can be an interesting laboratory for Luttinger liquid 
physics. These 1-1-1 dimensional systems have been extensi vely studied in the presence of correlated and point 
disorder 0, H, S [l^l ) or an isolated straight columnar pin [ll|, [12, [ll] ■ 

The continuum elastic energy for the coarse-grained displacement field u{x,t) for such a system of vortex lines in 
the presence of a straight columnar pin at x = and a transverse magnetic field H±^ is 
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where cn and C44 are the compressional and tilt moduli, no is the average vortex density which depends linearly on 
the external magnetic field, 0o is the quantum flux and is the pinning strength. 
A Luttinger liquid parameter g 
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can be defined from this long distance free energy. The long distance probability density derived from ([T]) exhibits 
Friedel-like oscillations, modulated by an exponential: 



{n{x)) - no 



cos(27mo2;) 



exp(-|x|/^). 



(3) 



The exponent a assumes different values for values of g above and below unity: a = 2g — 1 for g < 1 and a = g for 
g > 1. Here, ^ is a coherence length inversely proportional to the transverse magnetic field, ^ ~ H± or equivalently 
the relative tilt between the columnar pin and the direction H of the in-plane magnetic field of the slab. 

When the vortex array is dilute enough, or in the presence of a short range hard core repulsive interaction, the physics 
maps onto the free fermion problem. For a discussion in the context of commensurate-incommensurate transitions in 
adsorbed monolayers, see Refs [13, El, EE 113 ■ that case the compressional and tilt moduli are independent of the 
interaction potential details and assume values such that g — 1. The probability density distribution can be easily 
evaluated in this regime by use of the left and right eigenstates of the single particle non-Hermitian Hamiltonian 
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FIG. 1: Vortices confined to a planar superconductor. A linear defect (depicted with a thick grey line) that changes orientation 
in the superconducting slab acts as a pinning potential. The external magnetic field H is along the r dimension, which 
corresponds to _ffx = 0. 

where we have set fc^ = 1, T is the temperature and 7 is proportional to the line tension of the vortices. The 
simplest situation arises when the relative tilt between the defect and H remains constant. This condition applies 
to the straight defect case, i.e. a simple columnar pin, in which case ft.(r) cx 77 is independent of r. However, 
nanolithographic techniques could allow for the controlled fabrication of meandering linear defects, which motivates 
us to study the effect of a sudden change of direction or an abrupt termination of a defect trajectory within the 
sample. Because the one-dimensional Luttinger liquid is a critical system (with g-dependent exponents), there can be 
a striking response to such perturbations. 

The case of a single flux line interacting with a meandering linear defect has been studied in Ref. p^ . In this paper, 
we study the more experimentally relevant case of many flux lines interacting with a single kinked or terminating 
defect, confined to a thin superconducting slab. In Section |lT] we derive an analytic expression for the probability 
density distribution and imaginary current of flux lines in the presence of a bent pinning defect. These analytic 
expressions are used in Section IlIII to understand how the bent line defect perturbs the (l+l)-dimensional vortex 
configurations. In Sec. IIV[ we discuss potential experiments and defect lines that terminate. 



II. PROBABILITY DISTRIBUTION OF FLUCTUATING VORTEX FILAMENTS 

In a spirit similar to Ref. p^ . where the interaction of a single flux line with a meandering linear defect with 
trajectory Xo (t) was discussed, in this Section we derive an expression for the probability density of an array of flux lines 
interacting with an attractive delta function potential whose direction changes suddenly along the superconducting 
slab. To allow for a simple analytical treatment, we consider only piece-wise constant defect trajectories, and focus 
in particular on a defect consisting of two straight segments joined at an angle. Unlike the single vortex system, 
in the many flux line case there is no critical tilt that defines the critical point of a delocalization transition. The 
large number of extended states that are occupied even in the ground state in the many particle system wash out the 
delocalization transition of the single bound state. 

The short range attractive interaction between a vortex at position Xi and the pinning potential is approximated 
by V[xi{T) — Xo{t)] — —VQ6[xi — a::o(r)], where xq{t) is the pinning defect trajectory. The vortex-vortex interaction 
is assumed to be U[xi — Xj\ = Uo6{xi — Xj), with Uq 00. Upon assuming TV flux lines of length 2Lr confined to a 
length Lx with periodic boundary conditions at and (see Fig. [1} , we denote the positions of the vortices in the 
initial time-like slice as {x} = (xi, X2, ■ ■ ■ , xjy) and in the final slice at rj as {x'} — (x'l, x'2, ■ ■ ■ , x'^). The (classical) 
partition function for N flux lines then reads: 

Z = Z[{x'},Tf;{x},Ti;xa{T)] 

= / nPx.(r)exp(-i^ I drg (^^) -f ^[^^.(r) - xo(r)] -f ^ [/[x.(r) - x,(r)]] ) 
•'1=1 i i ^ i 



where 7 is the coarse-grained line tension. We continue to set Boltzmann's constant to fc^ = 1. 
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The transformation Xi(r) — > yi(r) = Xi{T) — Xo{t), i = 1, . . . , N for every vortex trajectory, enables us to switch to 
a frame of reference with a straight defect but varying externally imposed transverse magnetic field H± oc /i(t). We 
can now write the partition function of the many particle system as: 

Z = Z[{x' - Xo{Tf)}, Tf, {x - xo{T^)}, n; /i(t)] (5) 

where /i(t) = dxo{T)/dT, and the /i-dependent analogue of Eq. ([5]) is given by Eq. ^ 



Z^l f[Vx,ir) exp(-l ^ | dr [| (^^^ + ^r)) ' + V[x.(r)] + ^ U[x,{r) x,{r)]\ ). 



(6) 



As discussed in Sec. HI we consider the g = 1 Luttinger liquid limit and treat the vortex filaments as hard core 
bosons, or, equivalently (by means of a Wigner- Jordan transformation) as non-interacting fermions. 2J| However, we 
believe that the results will be qualitatively correct in cases where the average vortex separation a is large compared 
to the London penetration depth A which determines the range of the inter- vortex interaction potential U (x) , where 
X is the vortex separation. The problem of calculating the classical partition function ^ can be reformulated and 
recast in quantum language by mapping the thermally fluctuating vortices onto fermions propagating in imaginary 
time. In this mapping T h and 7 — + m 1]. The partition function then simply becomes the propagator 

^= ({x'-xo(r;)}|T.{e-/- ^-^M/^}|{2;-2;o(r,)}) (7) 

where TV is the time ordering operator and Ti-ir) is the time-dependent Hamiltonian of the many particle system 
associated with Eq. ([6]) above 

nir) ^~Jdx (^^^L (J^^ i'R + r/i(T)^i^^fl,^ - K)^L(0)^fl.(0). (8) 

This Hamiltonian is non-Hermitian, the non-Hermitian term being generated by the galilean transformation in imag- 
inary time. We find it therefore advantageous to work with a non-Hermitian set of particle field operators 

i^R{x) = ^0^(x,/i)4. (9) 

k 

(f>L{x, h) and (t)\{x^ h) are the left and right, respectively, eigenstates of the single particle Hamiltonian (i.e. neglecting 
interactions) for h =const indexed by a wavevector k [19, 20j. Note that, with this definition some usual relations for 
the particle field operators do not apply, as ijj^j^ ^ tpji. However, anticommutation relations, such as: 

{Mx),Mx')}^6{x-x') (10) 

do hold, and the vortex density is given by 

n{x) = ^]l{x)^Pr(x). (11) 

We will chose periodic boundary conditions in the x dimension for (j>'l{x, h) and 0|,(a;, ft.), a choice which is appro- 
priate when N is odd. [25'] In the absence of a pin, the density of vortices is a constant {n{x)) = no = N/L^- When a 
(possibly T-dependent) pinning potential is present, the probability density reads: 

where [^'(t))^ is the result of the evolution operator acting on the initial condition: 

|*(r));j-r,{e-^-V«(-')d-'}|vI/^)^, (13) 

and similarly for l{^{t)\: 

^(^(t)! =i(vI//|T,{e-/"^«(-')d-'}. (14) 
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To simplify notation we now also set T = 1. Here, 1^'^)^ and are non-hermitian generalizations of the usual 

filled Fermi sea ground state in one dimension (see below). As in Fig. [2 is the half length of the slab in the 
time-like direction. 

To make the calculation of {n{x))r analytically tractable, we now consider the simple example of a defect trajectory 
with a kink at r = 0: /i(t) = /i8(— r) + /i'8(t), with two slopes h and h', 

xo{t) = hrei-r) + /i're(r). (15) 

Moreover, we assume ground state boundary conditions, namely: I^E"*) = \Gh)B. and {'^^ — L{Gh'\, where \Gh)R is 
the right N particle ground state of ^{t) for /i(r) = h and L{Gh'\ is the left N particle ground state of Hir) for 
/i(t) = h' . The above assumption is not unreasonable if we take L,- ± r 3> 1 so the system has "time" to relax to the 
ground state both before and after it approaches the kink. Here \Gfi)R is the filled fermi sea ground state constructed 
from the eigenvalues of Ref. [l^. It includes the bound state (see Ref. [19,]), represented as a dot on the Ree axis of 
the energy spectrum diagram in Fig. HUja) and the extended state that occupy the paraboloid part of the spectrum. 
To compute the probability density distribution of the flux lines for r > 0, we rewrite the time evolution operator 

as 

Tr'{e-^--^ -i-' n-')} ^ W(r'>0)g-L. •H(r'<0) (^g) 

and insert a complete set of (normalized) many body eigenstates / = \i^h')R L{Kh'\- After some straightforward 
algebra we get: 

where 7?, is a normalization constant equal to 

n = L{Gh'\Gh)R, (18) 

and the sum is over all possible TV particle eigenstates of 7Y(r > 0). Here Ex{h') is the energy of the many body 
state \Kh')R, and similarly Ecih) is the energy of L{Gh'\- Note that in general Exih) has both a real and imaginary 
part, since the Hamiltonian it corresponds to is non-Hermitian and that although (|-R';i)_r)^ L{Kh\, these left and 
right eigenstates have the same eigenenergy. 

Since h{x) is a single particle operator, the matrix element L{Gh'\n{x)\Kh') r is zero for excited states with two or 
more particle-hole excitations. The summation is over all Khi states, where 



27r 27r 27r 

\Kh')R = I - kf, -kp + — , kh - — , fcp, kh + — , kp). (19) 

Ijrf. Ljrf Ijrf 



kp is the Fermi momentum, which equals kp = -f-iN — 2) ~ tthq for h < he- Although not explicitly shown in this 
notation, the bound single particle eigenstate is also included in the particle-hole excitation summation. These states 
correspond to single particle- hole excitations of the filled fermi sea \Gh')R. as shown in Fig. [Til 

Eq. (|17p effectively becomes a summation over all hole momenta \kh\ < kp and all particle momenta |fcp| > kp. 

The non-Hermitian particle field operators defined earlier lead to the matrix element: 

LK<^h'\nyx)\Kh')R - I for Kh' - Gh' 

where 4>^{x, h'), 4>\{x, h') are the single particle wavefunctions of the particle-hole pair and the n summation is over 
all occupied single particle eigenstates in Gh- It is easy to see that L{Kh'\Gh)R = det(C'), where the matrix C has 

elements Ci fc = / dxcpf {x, h')(j)^{x, h). 

L{Kh'\Gh) R. becomes negligible for large kp momenta, which allows us to introduce a cutoff kc at the summation in 
Eq. ^T7\ . An example is shown in Fig. [31 where we plot the logarithm of L{Kh'\Gh) r for a non-Hermitian field that 
changes from h — 0.05 to h' = 0.5 and an average vortex density Uq = 0.25. The value h = 0.05 was chosen instead 
of /i = for the r < tilt, due to the implicit assumptions used in the evaluation of the single particle eigenstates 
that appear in Eq. ([9]): To allow for analytical results in closed form, (f)\ and </)^ were derived under the assumption 
^-L^hj/T ^ 2^ which for finite systems is violated as /i ^ 0. The small initial tilt does not qualitatively change our 
observations. 
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(b) 

FIG. 2: Energy spectra and momenta for many body states for N oo when h < he = Vo/T. The isolated dot on the Ree axis 
represents the bound single particle state and does not appear on the wavevector line, but can participate in a particle-hole 
pair, (a) Ground state, (b) Single particle-hole excitation state. 
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FIG. 3: Logarithm (base ten) of the absolute value of the matrix element L{Kh> \Gh)R where K^/ differs from the ground state 
Gill by a single particle-hole excitation with hole momentum < kp and particle momentum |A;p| > kp- L{Kiii\Gh) r is 
maximum for \kp\ — » fcj and at kp ~ Skp it has decreased by more than 2 orders of magnitude. 



Having calculated {n{x))r, we can undo our transformation, x ^ x — xq{t) to get the actual probability density in 
the original frame of reference: 

{n{x))r = ^ 0L(^-a^o(T),/i')'^^(^-a;o(T),/i')(Lifft'|G/.)/?e-"(^-('''^-^'^("'»+n„,(x-xo(T)) (20) 

kF<\kp\<kc 
\kh\<kF 

tp{h') — e/i(/i') is the particle-hole energy and n/i'(x) is the vortex density {n[x,T)) for constant /i(t) = h' . 
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h=0.02 ^ h=0.4 
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FIG. 4: Density of vortices when the tilt of the hnear defect changes from h =0.02 to h' =0.4. For better visuahzation, the 
range of the n{x,T) values was restricted to [0.16, 0.27]. On the right: snapshots for times r =-30, r =0 and r =48. Note the 
Friedel oscillations. The dashed curves represent the exp(— |a;|/^)/|x| envelope of Eq. ([Sj. Close to r = 0, an assymetry develops 
in the n{x,T) profile, and therefore such an envelope cannot be defined. Lengths are measured in units of [x] = T'^/{'yVo) and 
imaginary time in units of [r] = T'^ / {jVq) . 



We can easily adapt the above treatment to get the probability density for r < 0: 



kp<\kp\<k^ 
\kh\<kF 



Another quantity which we can give us some insight to the properties of the transition between the h and h' sections 
of the tilted defect is the expectation value of the current operator [l^, [ill : 



dH{T) _ iT 
' dh ~ T 



dx 



(22) 



Note that although the Hamiltonian depends on h, the current operator itself is explicitly /i-independent. With our 
definition, J is not a current density operator, but an integrated, position independent quantity, a measure of the 
total transverse magnetization due to the tilted flux lines. 

In the absence of a pinning defect, the current (J) is zero since the flux lines are parallel to the inducing magnetic 
field. When a defect is present, the flux lines bend in the vicinity of the tilted defect, thus creating a non-zero 
/i-dependent current. As long as the magnetic field remains parallel to the r-axis, the total current (J) is independent 
of the system size, since vortices far from the defect are parallel to the magnetic field and do not contribute to the 
current. 



III. OSCILLATIONS IN THE TIME-LIKE DIMENSION 



In this Section we present some results obtained using the formalism discussed in Sec. |TT1 and examine the flux line 
conflguration in the vicinity of the defect kink. 

Figures m and [5] show the vortex density distribution for a system of = 91 vortices penetrating a planar slab of 
length = 360 (measured in units of T'^/jVo), and interacting with a linear defect with tilt h = 0.02 for r < and 
h — 0.4 and 0.9 respectively for r > 0. 

As T — ^ cxD, the probability distribution is peaked on the defect and exhibits oscillatory behavior with amplitude 
that decays as a power law modulated by an exponential, as described by Eq. ([3]) with a = 1 ("Friedel" oscillation 
behavior [ll|, [13 )• For h ~ the vortices form a periodic array, centered on the defect. The situation changes 
drastically close to the defect kink, especially for r > 0. In this case there is a competition between the externally 
imposed magnetic field H, which alone would result in an array of vortices parallel to the field, and the tilted defect. 
The vortex which was localized on the defect for r < follows the defect until forced out by another vortex, which 
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h=0.02 ^ h=0.9 
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FIG. 5: Above: Density of vortices when the tilt of the linear defect changes from h =0.02 to h' =0.9. For better visualization, 
the range of the n(a;,r) values was restricted to [0.16, 0.27]. Below: snapshots for times r =-30, r =0 and r =48. The dashed 
curves on the top and bottom graphs represent the exp( — |a;|/^)/|a;| envelope of Eq. Q. ^ is inversely proportional to h, so for 
high tilts no Friedel oscillations can be seen. Lengths are measured in units of [x\ = r^/(7Vb) and imaginary time in units of 
[r] = T'/i^Vi). 




FIG. 6: A probable configuration of flux lines on the planar superconductor. In this sketch, flux lines have been drawn to 
roughly follow the probability density maxima. The dotted ellipses demonstrate the likely position on the defect for the vortex 
exchange. 



takes its place. Since the average vortex- vortex distance is a = 1/no, the localized vortex stays on the defect over a 
time- like distance At ^ l/(no/i). The exchange of vortices that are localized on the defect (one vortex enters while 
the next escapes the defect) takes place at the neighborhood of the probability density local minima n(a;o(T(), rt). 
However, as r — > oo, n(a;o(Tt), rt) becomes constant so the coordinate tj of the exchange is well defined only close to 
the kink. We expect that the r-independent probability distribution nh'{x, r — > oo) is approached exponentially fast. 
A typical vortex configuration is presented in Fig. [6l 
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FIG. 7: Total imaginary current (see Eq. (|22|l 'l for t > transition h — 0.05 — > /i' = 0.5. The dotted line is a plot of (n(a;o(r)))T-, 
the probability density on the defect in arbitrary units. Inset: Imaginary current for r < 0. 




FIG. 8; Logarithm of the amplitude describing the oscillatory decay of n{x, r) away from a bent defect for a transition 
h = 0.05 ^ h' = 0.5, with = 91 vortices and = 360. The length and time units are [x] = r^/(7Vb) and [r] = T^/{'yVo ) 
respectively. The frame of reference has been shifted so that the defect (and the amplitude maximum) remains at a; = for 
both positive and negative t. 

The kink breaks time-translational invariance and leads an oscillatory behavior in the current, as well as the density, 
with period At' = l/{noh'). The imaginary component of the "current" is maximized when the vortex follows the 
defect, and minimized at T(, the exchange point. The oscillations in the current die out far from the kink, as the 
"jumps" now occur with equal probability at all times t and a traffic jam profile similar to the one described in 
Refs. (ill , [l^ is formed. As can be seen in Fig. [3 the maxima in the current coincide with the maxima of n(xo{T), t). 
An oscillatory behavior of the current also appears for r < as shown in the inset of Fig. [3 this time with period 
At = At diverges as ft. — > and no oscillations in imaginary time are present. 

The length of the superconducting slab necessary to observe the phenomena discussed in this work would primarily 
depend on the imaginary time vortex-vortex collision period Ar' ~ a/h' . If we take U{x) to be proportional to the 
modified Bessel function of zeroth order (see for example and references therein), then for a ~ 3.5A we have 
U {a) /U {X) < 0.05, so we would expect that the behaviour of the vortex system at that density could reasonably 
be approximated with a non-interacting model. For A ^ 40nm, then for a system of 10 vortices we would need a 
superconducting wafer of width ~ 1.5/im and length ~ 1.5/im//i', well within the experimentally accessible region. 

The behaviors described before are summarized in Fig. [8l where we plot the logarithm of a local amplitude A(x, r) 
of the oscillations for a system of = 91 vortices (again spread over a length L^. = 360), in the presence of a defect 
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Number of vortices 
on concave side of defect 
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FIG. 9: Total number of vortices iV^ = J^^^^J^'^^''^^ {n{x))^dx on the concave side of the defect for a transition h = 0.02 — » 
h' = 0.7. 

with h = 0.05 for T < and h = 0.5 for t > 0. The amplitude is a coarse-grained oscillation envelope obtained by 
subtracting the density of a local minimum from the density at the previous local maximum: A{{x^in + a;max)/2, t) = 
n{xraa,XT t) ~ '^■(a^min, t) . For casicr visualization, in this picture we have shifted the frame of reference to one which 
moves with the defect, so that the center of the defect is at a; = 0. As always, the defect changes slope at r = 0. 

The picture consists of an almost uniform in x oscillatory background with a maximum at r = which decreases 
exponentially fast in r (appearing as a straight line in a logarithmic plot) far from the kink. [2^ The oscillatory 
background is solely due to the presence of the kink in the defect trajectory. This abrupt change results in enhanced 
positional order of the vortices far from the defect. The "tent" -like structure emerging from the background as 
r ±oo is due to the oscillatory behavior of the time independent part of the density nh{x) for r < and rih' for 
T < 0. As the jump from /i to /i' at t = is approached, the width of the "tent" -like structure, which is a measure of 
the coherence length f ^ decreases and the "tent" becomes narrower. Of course we should bear in mind that the 
quantity A(x, r) is not well defined close to x — 0, where n{x, r) changes much faster than A ^ L^/N, the wavelength 
of he Friedel oscillations. With this caveat, the oscillatory behavior at a; = with period At for r < and At' for 
T > can be understood as a signature of the vortex interactions mediated by the kink in the defect. 

Far from the kink, at each time slice the vortex density distribution is symmetric about the defect position xq{t). 
Although not evident at the resolution of Fig. [51 the reflection symmetry of n{x, r) across the defect is broken near the 
kink is broken. The density of vortices is higher on the concave side of the defect, although only by a small amount, 
as seen in Fig. [9] This effect was also present in the single vortex- meandering defect system, and disappears when 
h' ^h^. 

The average vortex positions {x)i (corresponding to the local maxima of the vortex density) far from the defect 
kink follow the well known Friedel oscillation pattern [ll|, forming a periodic array with lattice spacing I/uq. 
The phase of the oscillations far from the defect, depends on the tilt h and the strength of the attractive potential. 
However, as the defect kink is approached the simple periodic pattern in the vortex positions breaks down, and a 
phase shift is introduced. To explore how the defect kink affects {x)i in Fig. [TUTa) we plot the local maxima (average 
position of each flux line) of the probability distribution at imaginary time t = (squares) and at t = 60 (diamonds). 
Each column corresponds to a different value of h' . In all cases, h = 0.02. The blue cross corresponds to the position 
of the defect at t = and the red cross at t = 60. 

At each column, the data points beetween two crosses correspond to the vortices whose trajectories crossed the 
defect between t = and t — 60. By simple inspection it can be determined that there is a minimal change in the 
average vortex position between t = and t = 60 for vortices that did not cross the defect. However, vortices that 
cross the defect get trapped and move with it for some time until forced out by the next vortex that impinges upon 
the defect. 

A more detailed description is given by Fig. llOf b). where we plot the shift in average position of the vortices {x)i 
versus vortex number i for different tilts. We arbitrarily labeled the vortex at distance (x) ~ 50 as vortex 1. A positive 
A{x)i corresponds to a vortex position shifting to the right at t = 60 with respect to its original position at t = 0. 
For all tilts h', the vortex which is immediately to the left of the defect is shifted to the left. As discussed previously, 
vortices enter, become trapped and then exit the defect at periodic intervals. For certain tilts h' , the imaginary time 
slice T = 60 happens to be when a vortex enters the defect, by abruptly turning to the left. This is the origin of the 
maxima that can be observed on Fig. [TUf b). 

Ignoring those t dependent features, all displacements A{x)i of the vortices that cross the defect collapse on the 
same value, a trapping displacement /tr approximately equal to Ztr — 0.38a, where a is the mean vortex spacing. The 
trapping displacement does not depend on the tilt h' , as similarly defined quantities for single vortex systems, i.e. the 
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FIG. 10: (a) Average position flux lines at r = (squares) and at r = 60 (diamonds). Each column corresponds to a different 
value of h' . In all cases, h — 0.02. The blue cross corresponds to the position of the defect at r = and the red cross at 
r = 60. (b)Shift in average position of the vortices {x)i versus vortex number i. The vortices that cross the defect exhibit an 
approximately constant shift in position, idependent of h' . 



trapping length in [19] and the vortex shift m{h) in |18| . do. This contradicts the intuitive expectation, derived from 
the single pinned flux lines that ^tr should decrease as h' increases. 



IV. POTENTIAL EXPERIMENTS AND DEFECT LINES THAT TERMINATE 



The above results and discussion about the probability density are applicable only far from the edges of the slab, 
a region difficult to probe experimentally. The use of ground state initial and final boundary conditions to derive 
the probability density in the core of the slab requires that the coordinate r of interest be far from the boundaries. 
The experimentally relevant boundary condition (final or initial state) is an extension of the single particle boundary 
conditions f dx(a;| for many particles, namely (^'■'^1 = / dxi J^^ dx2 ■ ■ . J^"""^ dxN{xi, X2, ■ ■ ■ ,xn\ (and similarly for 
the initial state). If we examine the probability density distribution at the upper boundary (which describes the exit 
positions of the vortices ), we can assume ground state boundary conditions at the other boundary provided L^- ^ 1- 
The energy spectrum of the single particle states has an imaginary component. Due to this component, irrespective 
of the exact form of the exiting boundary conditions, we expect to observe oscillatory behavior with respect to the 
distance between the kink and the exit surface Lj. The probability density at the upper boundary reads: 
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FIG. 11: Sketch of possible experimental setup for the kinked defect system. 
The MFM tip probes the neighborhood of the defect exit position. 



The defects are represented with dark gray lines. 



where AE{Kh') = E{Kh') — EiGw)- If, however, L/ ^ 1, then we can approximate the sum in the denominator of 
Eq. ([221) by the term Kh' = Gf and find for the density of exit points: 



{n{x))Lf ^ n^j{x) + 



1 



{^f\Gh'){Gh'\Gh) 



J2 i{^f\Hx)\Kh')-n^fix){^f\Kh'))e-^''^''-'^'^f{Kh'\Gu) (24) 



Here, nq,f{x) is the probabiUty density distribution at the boundary for Lf — > oo: 

_ {'ff\h{x)\Gh') 



n^f{x) 



{■^^Gh 



(25) 



For sufRciently large Lf, the sum in the second term of Eq. (|24|) is dominated by the lowest energy eigenstates Kh'. 
These are the single particle excitations with energy ReE{Kh/) ~ N/L^. and momentum Ap ~ 2kp (the equality 



holds for La, oo). Note that e"^-^ ^f'^ ~ ^-i2iikFLj ^ ^-z2^nohLf ^ r^j^^^^ j^^^p Ax = x - xo{t), the distance 
from the defect, constant we expect to observe a periodic modulation in the r direction of the vortex probability 
distribution. This periodicity should manifest itself in magnetic force microscope (MFM) experiments as it traverses 
the slab [13, [1^ , since the force necessary to pull the vortex from its exit position would depend on the length L/ of 
the tilted segment of the defect. For example many linear defects with varying tilted segment length can be etched 
on the same superconducting slab, as in Fig. [TlJ If the defects are sufficiently far apart so that we can assume that 
each vortex interacts with only one defect, then we expect 7i(x,t) of the single defect system to approximate well 
the density at the neighborhood of each defect. With the proper choice of the tilted segments of the linear defects 
imprinted on the slab, one should be able to observe a sin(27mo^L/) dependence of the probability density at a fixed 
distance from each defect end point. 

An interesting question that can be explored through a similar approach is the case of a defect terminating while 
still inside the slab. In this case the localised bound state that existed for h smaller than he, the critical tilt for the 
delocalised transition disappears for r > and the spectrum consists only of extended states. The localised vortex 
starts diffusing as it approaches the defect end and spatial phase information is quickly lost fFig. [T^. The "locking" 
of phase in the time-like dimension that was observed for kinked defects and the oscillations in r of the vortex density 
that this resulted to are not present in the case of the terminating defect. This can be seen in Fig. [T31 where the 
defect is at tilt h = 0.7 with respect to the externally imposed field. Vortices still enter and exit the defect forming 
the traffic jam discussed in Refs. [ll[ and [12], but the exit position of the last trapped vortex can vary so oscillations 
of the vortex density in the time like dimension are not observed. 
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APPENDIX A: DERIVATION OF THE NON-HERMITIAN HAMILTONIAN 



The mapping of the classical flux line system to imaginary bosons (or fermions) is explained in many publications, 
see for example Refs. [isl [loj . However for completeness, we briefiy review it here. The partition function for one 
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FIG. 12: Density of vortices with a terminating linear defect with tilt h =0.02. The defect terminates at r = 0. On the right; 
snapshots for times r =0, r =2 and r =6. The slight assymetry of the density profile near r = is due to the nonzero tilt h of 
the defect, and disappears for h ^ 0. The density distribution delocalises exponentially fast after t = 0. Lengths are measured 
in units of [x] = T^/{'yVo) and imaginary time in units of [r] — T"^ /(jVq). 
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FIG. 13: Density of vortices with a terminating linear defect with tilt h = 0.7. No oscillatory behaviour in imaginary time is 
present. 



flux line interacting with a meandering defect is: 



Z[xf,Tf;x,,n;h{T)]^ / I?x(T)e-^[-Ml/^ 

J x(Ti)—Xi 



where 



E[x{t)] = J dr 



7 / dx[T) 
2 V dr 



Vo 6{x - xo{t)) 



(Al) 



(A2) 



Xi and Xf are the positions of the vortex at T,j and Tf respectively. A change of variables yir) = x{t) — Xq{t) will 
transform Eq. (|Aip to a functional integral over the variable 2/(t) with energy: 



E[y{T)] = j dr 



7 ( dyjr) 
2 V dr 



+ h{T)] -VoSiy) 



(A3) 



where hir) — dxa{T)/dT. 

The partition fmiction obeys the equation: 



d T'^ d"^ d 

-Q:^2[y,T;y„Ti;h{T)] = ( - - '*('^)^^ ~ VQS{y))Z[y,T;y„n;h{T) 



(A4) 



Thus, Z[x, t; yi, Ti] h{T)] can be thought of as the quantum mechanical propagator of a particle propagating in 
imaginary time, with non-Hcrmitian Hamiltonian: 

^(^) = ^P'' - Mr)p ~ Vo6{x) (A5) 
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where 73=?^. The generahzation to many Hnes is straightforward. 

The non-Hermitian single particle eigenstates 4>fi{x, h) used throughout this work are derived from solving Eq. (|A5|) 
for /i(t) = cons[l^. When h < he = Vq/T, for a system of length Lx, with x G [0,iva;] and periodic boundary 
conditions, the (unnormalized) extended states read: 



Q-~ikx-\-6Kx ^ikx — 2'yhx/T — SKx 

-ikLj,-^5KL^ _ 1 pikL^—2'yhL^/T-5KL^ 



(A6) 



where 



6k = 



k + i^{h- 



T I 



(A7) 



provided that e 



<C 1. The bound state similalry reads 



j_(Yo. 



g ^ T" h^X 



(A8) 
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